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Ahatnet 

In  this  paper  we  present  a  general  model  of  the  completion  time  of  a  single  job  on  a  computer  sys¬ 
tem  whose  state  changes  according  to  a  semi-Markov  process.  When  the  state  of  the  system  changes  the 
job  service  is  preempted.  The  job  service  is  then  resumed  or  restarted  (with  or  without  resampling)  in  the 
new  state  at,  possibly,  a  di^frent  service  rate.  Different  types  of  preemption  disciplines  are  allowed  in  the 
model.  Successive  aggregation  and  transform  techniques  are  uaed  to  obtain  the  Laplace  Stieltjes 
Transform  of  the  job  completion  time.  We  specialise  to  the  case  of  Markovian  state  process.  Finally,  we 
demonstrate  the  use  of  the  techniques  developed  here  by  means  of  an  ^>plication. 


*  Thb  work  wm  npporud  is  part  kjr  ibc  Air  Foret  Oniet  of  SeitBUfie  Rtsesrek,  hj  the  Armj  Rtstarch  OfTiet  isdcr  costraet 
DAAG3S-S44I045  sad  b;  Ibt  Natioa^  Sdtaet  Poaadatioa  sadtr  graat  MCS-UOSOO. 
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1.  Imlnimetimi 


A  computer  system  may  change  its  state  (e.g.,  its  mode  c^raUon)  due  to  different  events.  The 
system  behaviour  can  adequately  be  described  by  a  stochastic  process,  referred  to  as  the  "structure-state 
process*.  Consider  a  job  to  be  processed  on  such  a  system.  The  job  service  is  preempted  due  to  changes 
in  the  state  of  the  ^stem,  and  may  be  continued  ot  repeated  at,  possibly,  a  different  service  rate. 
Clearly,  the  job  completion  time  (i.e.,  the  total  time  spent  in  the  qrstem  to  complete  the  job)  is  affected 
by  changes  in  the  system  operation.  It  is  important  to  distinguish  different  types  of  service-preemption 
interaction.  In  the  preemptive-resume  (prs  )  discipline,  the  job  service  is  resumed  from  the  point  where  it 
was  preempted.  In  the  preemptive-repeat-identical  (prt )  discipline,  the  job  is  restarted  from  its  beginning. 
In  the  preemptive-repeat-different  (prd)  discipline,  the  job  is  restarted  with  a  new  work  requirement 
which  is  statistically  independent,  and  identically  distributed  to  the  original  work  requirement. 

Several  authors  have  stpdied  the  job  completion  time  in  special  cases.  Caver  [4]  and  Nicola  [9]  con¬ 
sidered  a  single  server  system  subject  to  different  Qrpes  of  Poisaoii  interruptions.  In  their  system  a  job  is 
serviced  only  in  one  state  at  a  constant  service  rate.  Castillo  and  Siewiorek  [2]  presented  a  model  with 
two  types  of  Poisson  breakdowns  and  repair.  Puri  [10]  studied  the  distribution  of  the  cumulative  service 
in  the  case  where  all  preemptions  are  of  the  pr$  type;  Uiis  is  shown  to  be  dual  to  the  distribution  of  the 
completion  time  of  a  given  job  [7].  Under  Markovian  assumptions  concerning  the  system  changes  (i.e.,  a 
Markovian  structure-state  process),  we  presented  a  detuled  analysis  in  the  cases  where  different,  pure  and 
mixed,  types  of  preemptions  are  allowed  in  the  model  (6,7).  In  [5],  Iyer  et  al.  considered  the  computa¬ 
tional  aspects  and  developed  a  procedure  for  the  computation  of  the  moments  of  the  cumulative  service  in 
the  case  where  all  preemptions  are  of  the  prt  type  and  the  structure-state  process  is  Markovian. 

In  this  paper,  we  extend  the  results  presented  earlier  to  allow  semi-Markovian  system  changes. 
Furthermore,  all  types  of  preemptions  may  be  present  in  the  same  model.  This  is  a  general  framework 
that  includes  all  previous  models  as  special  cases,  and  extends  the  results  to  permit  the  analyns  of  fairly 
complex  systems.  In  section  2,  we  describe  the  mathematical  model.  The  job  completion  time  is  analysed 
in  section  S.  We  specialise  to  the  case  yrhere  the  structure-state  process  is  Markovian  in  section  4.  In  sec¬ 
tion  5,  we  demonstrate  the  use  of  the  techniques  developed  here  by  means  of  an  application. 


2.  TUBa^UUd 


Crainder  »  job  with  work  reqairement  ,  to  be  processed  on  a  computer  system.  The  work  require¬ 
ment  is  measured  in  work  units,  e.g.,  number  of  instructions  to  be  executed.  Assume  that  Bum  positive 
random  variable  with  cumulative  distribution  function  G(g)^  P(B  <  *).  The  computer  system  can 
be  described  by  its  structure-state  process  {Z(t),t  >  0),  which  is  assumed  to  be  a  continuous  time  sto¬ 
chastic  process  deHned  on  the  state  space  {0,1,2,. ..,a  }.  The  woric  rate  of  the  computer  system  is  r,-  >  0 
(units  of  work  per  unit  time)  when  it  is  in  state  i.  The  state  0  is  an  absorbing  “failure”  state  with 

ro  0.  Each  state  in  {1,2 . n  }  is  assumed  to  be  prs  or  pri  or  prd ,  as  defined  in  the  previous  section. 

Let  St(Sg,St)  be  the  set  of  all  prs  (pri.prd)  states.  Thus  Si,S2,St  is  a  partition  of  {l,2,...,n}.  We 
assume  that  the  structure-state  process  is  stochastically  independent  of  the  work  requirement  of  the  job. 

Every  time  the  system  virits  a  prt  or  prd  state  all  the  previous  work  b  lost  and  the  job  is  restarted. 
Define  T  (x )  to  be  the  amount  of  time  needed  to  complete  a  job  with  work  requirement  of  x  units;  and 
T  to  be  the  time  needed  to  complete  a  job  with  random  work  requirement,  B .  We  are  interested  in  the 
distribution  of  T  under  the  assumption  that  the  job  starts  being  processed  immediately  after  the  transi- 
rion  of  the  structure-state  process  to  the  initial  state. 


Dellne  the  distribution  functions: 


F.(t,x)  =  P(r(x)<  <  |Z(0)=.), 

F(<,x)  =  P(T(x)<  0. 

<  t  |Z(0)=.), 
F(t)~P{T  <  I), 


<  ,x  >  0,  1  <  »  <  « 
l,x  >  0 

1  <  «  <  * .  ‘  >  0 
I  >  0. 


Define  the  LST  s*: 


/;’(s,x)-E(e-*^‘')|Z(0)=i),  »  >0,l<i  <n, 

F  (s,x) *  >0 

i-i 


F/  (s )  -  1  Z(0)-.  )  -  fF/  (a  ,x)dG(x), 

0 


1  <  i  <  n 


('  )  StaotM  the  Lspiset  StieItjM  Trsaiforiii  (LST ). 


I'J  '•^.."l  I."!'.'*  ."»  ”11  '.'V:"ii  L  'i!"*! 


r  (,)^E(c-n^  E/’/(*)/‘(^(o)-0- 

<«i 

It  18  clear  that  Ff  (•  ,x )  and  Ff  (« )  are  the  key  quantities  in  the  analysis  of  T .  b  the  remaining 
analysis,  we  make  the  assumption  that  {^(t),t  >  0}  is  a  semi-MarkoT  process  (5A<P).  Let  be  the 
holding  (sojourn)  time  in  the  initial  state,  i.e. 

H  =min{t  >  0:Z (t )  Z (0)} . 

Let 


& 

i 


iU 

h 

w  • 

v' 


I 


e,y(x)^F{If  <  x-,Z[H+)=i  I  Z(0)=i},  0  <  ij  <  " . 

be  the  distribution  of  the  sojourn  time  in  state  t,  given  that  a  transition  to  state  j  took  place. 
Q{x)—  [<2,y  (z )]  is  called  the  kernel  of  the  SMP .  Let 


Q.(*)— EQo(*).  !<•<". 

i-o 

be  the  distribution  of  the  sojourn  time  in  state  • ,  and  deHne 

0 

and 

eo 

J  f^dQi(x). 

0 


3.  The  Anxigau  •/  Ihc  CsinffctMa  Ttmc 

In  this  section  we  describe  a  systematic  procedure  to  derive  the  distribution  F(( )  of  the  job  comple¬ 
tion  time  T .  The  method  of  analysis  is  best  described  as  "progressive  aggregation”.  To  begin  with,  we 
study  the  structure-state  process  restricted  to  move  only  in  5|,  then  we  study  it  restricted  to  Si  U  Sg, 
and  finally  we  study  the  general  case,  when  it  moves  freely  in  5]  U  U  Sy 


6 


3.1.  The  Analif$i$  of  the  Completion  Time  in  the  pro  Set  of  Stmtee 
Suppose  the  structure-state  process  is  initially  in  S|,  and  let 
t/i  —  min  {t  >  0:Z(t)^  Si). 

Thus,  Ui  represents  the  first  time  the  structure  state  process  visits  a  preemptive  repeat  state  or  the 
absorbing  failure  state.  Now  consider  a  job  with  w<»fc  requirement  s ;  its  job  completion  time  is  r(z ).  If 
^  ^  ^1  job  completes  while  the  structure-state  process  is  still  in  5],  if  ^(z )  >  [/|  the  job  has  to 

be  either  restarted  in  a  state  in  5^  U  5s  or  the  job  is  never  completed  due  to  an  absorbing  failure.  Define 

3/;.,(s,z)  =  £{e-"»'<«>;r(z)<  U,  |  Z(0)=.),  .  €5,, 

to  be  the  LST  of  the  job  completion  time  T(x)  when  the  job  completion  takes  place  before  the 
structure-state  process  leaves  the  set  5|.  Notice  that  ,z)  is  the  LST  of  a  possibly  defective  ran¬ 

dom  variable,  sinee  ^(^(z)  <  Ui  |  2(0)’^i)  <  1. 

Define,  for  (  C  Si  and  j  ^  5|, 

Ml],  Jo  ,z )  -  £(<’“'*;r(z)  >  Ui,Z(Ui+)^j  |  Z(0)  -  • ),  •  €  5„  /  ^  5,. 

When  the  structure-state  process  starting  in  t  €5|,  leaves  5|  before  the  job  completes  and  enters  state 
j  ^  5j,  the  LST  of  the  total  time  spent  in  the  set  5i  is  given  by  Afi(j(e  ,z )  .  Now  let* 

00 

M I'fKs  ,w )  =  /  e"**  Af  i.,(s  ,z  )dz  ,  i  €  5, 

0 

and 

00 

A/ (s  ,«,)  =  /«—  Mi',  Je  ,z  )dz ,  .•  G  S„  j  ^  S,. 

0 

The  theorems  below  give  the  equations  to  obtain  the  above  quantities. 

Theorem  1.  The  double  transforms  A/i,,*(«  ,w ),  i  G  5 1,  satisfy  the  following  equations 

~ E  9a  i  G  5,.  (1) 

•+^«’  *C5, 

Proof:  Conditioning  on  H ,  the  sojourn  time  in  the  initial  state,  we  get 


(•)  dtDotn  tht  Lspiset  trsatrorm. 


e  if  A  >  ar  /r,- 

ies^ 

UnconditioBing  <m  /f ,  we  get 


•  A,  *€Si  0 

•  A, 

=  e ■“  (l-«.  (x  /r, ))  +  E  /  e M,;*(e  .x -r.-  A  )</(?;*  (A ). 

t€5,  0 

Multiplying  both  sides  by  e'**  and  integrating,  we  get 


00  OO  •  /'l 


t€S,  0 


which  yields  equations  (1). 


Q.£.D. 


Tkeorem  t.  The  double  transforms  Af  ( ,w),  %  C5|,  y  ^  Si,  satisfy  the  following  equations 

^ i’u(«  >«»)**  -ir  <?.7  (•  +f<  «> )  +  E  (•  +'( »» ).  • 

*  *€S, 

Proof:  Conditioning  on  ff ,  the  sojourn  time  in  the  initial  state,  we  get 


E(c\T(z)  >  l/„Z(i;,+)-  ;•  |f/=A,Z(0)=.) 


fe  ,  ifA<x/r, 

~  1  e  ■**  Af  1,*  ,,(s  ,x  -r,-  A  ),  if  A  <  x  /r, 


ifA<x/r,-  andZ(A)=y^5i 
if  A  <  X  /r,  and  Z  (A  )=A  G  5  j 


Unconditioning  on  /f ,  we  get 


*/f,  »/', 


*€5,  0 


Multiplying  both  sides  by  e'”*  and  integrating,  yield  equations  (2). 


Q.E.D. 


Remark:  If  the  structure*8t8te  process  is  a  CTMC,  the  above  theorems  yield  propositions  5.1  and  5.2  in 


Equations  (1)  and  (2)  can  be  solved  to  get  Mij\t  ,w )  and  Af  i,,*y(s  ,w )  ,  i  €.5 1  ,  j  ^5 1  ,  which 
can  be  inverted  with  respect  to  w  to  obtain  Af,-j(s,x)  and  ,z),  respectively.  If  the  structure- 

state  process  is  a  CTMC ,  this  inversion  can  be  carried  out  rather  easily  [7]. 

The  quantities  ,x)  and  ,z)  completely  describe  the  job  completion  process  while  the 

structure-state  process  is  in  5  ].  Using  these  quantities  we  can  analyze  the  job  completion  process  while 
the  structure-state  process  is  in  5|2=Si  U  • 

3.2.  The  Analpeis  of  the  Completion  Time  in  the  prt-pri  Sett  of  Statet 

Suppose  the  structure-state  process  is  initially  in  5)2  (i.e.,  in  5]  U  ^a),  and  let 
Ux2  =  min  {t  >  O.Z (1)^5 la}- 

Thus,  U 12  represents  the  time  until  the  structure-state  process  visits  a  prd  state  cv  the  absorbing  failure 
state.  Now  consider  a  job  with  work  requirement  x ;  its  completion  time  is  T  (z ).  V  T{x)  <  Ui2  ,  the 
job  completes  while  the  structure-state  process  is  in  Sif.  If  T(x)  >  Un,  the  job  has  to  be  either  res¬ 
tarted  in  a  state  in  8%  with  a  resampled  work  requirement  or  the  job  is  never  completed  due  to  an  absorb¬ 
ing  failure.  Define 

Af,'2,  («  ,*)  -  £(«-*’■<' »;r(z)  <  U,2  I  Z(0)  =  f ),  .  €5,2  , 

to  be  the  LST  of  the  job  completion  time  T  (z )  when  the  job  completion  takes  place  before  the  process 
leaves  5,2.  For  •'  €5,2  and  /  ^5,2  ,  define 

•*;r(z )  >  t/,2,z(i/,2-(-)  =  y  U{0)  =  .),  .€5,2,  y^5,2. 

When  the  structure-state  process,  starting  in  .'  €  5 ,2,  leaves  5 ,2  before  the  job  completes,  and  enters  a 
state  •  ^  5,2,  the  LST  of  the  total  time  spent  in  the  set  5,2  is  given  by  M,2,i,y(e  ,z )  .  Now  let 

(s  ,z)  =  c'"^'-(l-(?,(z/r,))  +  E  /  .T-rih)dQ^(h  ),  .€52 

ACS,  0 

and 

A,;(s,z)=  /c--dg.y(M+  E  /  (A). 

0  0 


i,3  €52. 


The  next  theorem  provides  a  method  of  computing  the  LST  s  Mni(s  ,s ),  i  C5i2. 


Theorem  S. 

a)  The  conditional  LST s  ,x ),  t  GS3,  satisfy  the  following  equations 

Ml2.i {*.*)=  9/  (*.*)+  E  *0  (*  »*  Wli.A*  •*  )>  2-  I 

b)  The  conditional /.STs  A/i2_,(«  ,ar ),  1  C5,,  are  given  by 

A/ =  •(«,*)+  S  A/i,. •,/{«, s)3/,2.y(»,^),  i£5,.  , 

ycSj  ' 

Proof. 

a)  Let  I  C52  be  the  initial  state.  Conditioning  on  the  sojourn  time,  H ,  in  the  initial  state,  we  get 
£(«-*■<•)  ;r(x)  <  Uu\Z{0)=i,H=^h) 

c”**^'*,  if  h  >x/ri 

«'**  5^  A/i2,*(s  ,* ) 

*€5j 

+  «-**  2  (A/l’*(«,*-^A)+  E  A/i'*.y{s,*-r,  A)M,^,{#,x)),  if  A  <  x/r,-  . 

*eSj  yes. 

Unconditioning  on  H  and  rearranging  yields  equations  (3). 

b)  Let  I  €5]  be  the  initial  state.  Then 

£(e-*^<');r(x)<  C/,2U(0)=«) 

=  £(e-'n.);r(x)  <  U,  I  Z(0)  =  1)  +  £(e -*»•<* );t/,  <  TCx)  <  f/,2  |  ^(0)  =  t). 

It  follows  that 


M,2..(«  .*)  =  .*)  +  E  £(e''''‘;7’(*)  >  U„Z(U,+)  =  ;•  j  Z(0)  =  .  )  ,x) 

yes, 

which  yields  equations  (4).  Q.E.D. 

Remark:  If  the  structure-state  pTX>cess  is  a  CTMC ,  the  above  theorem  produces  theorems  5.1  and  5.2  in 


The  next  theorem  provides  a  method  of  computing  the  LST s  A/i2,,-,/(s  ,x),  i  G5i2,  j  First,  define 


•  /',  •  /'■ 

ai-  (s  ,x )  =  /  t "*  iQij  (A  )  +  E  /  « ■**  Mijt .,(s  ,* -r.  k  )dQit  (k  ),  i  G S,.  j  ^ 5 ,2 
0  *es^  0 

and 

»/r,  »A, 

(s  ,x )  =  f  dQn{k)  +  E  /  *  **  Af r^  )d(5,„  (A ),  ijk^S^. 

0  ■  €5]  0 

Tkeorem  4- 

a)  The  conditional  LSTs  A/,2,,  ,y(«  .x ),  t  GSj  .  j  ^ 5,2,  satisfy  the  following  equations 

A/i2,i.y(«  .*)  =  “tj(*  >*)  +  E  **(*  •*)’  *^^2.  y^5,2. 

b)  The  conditional  LSTs  A/,‘2,..,(*  ,i ),  1  GS,  ,  /  ^5,2  ,  are  given  by 

A/i2.<.,(<  ,3f )  =  A/,‘,-,,{s  ,1)  +  E  '*)*^»2.t.y(*  >*)>  •C5,,y^S,2. 

*€Sj 

Proof. 

a)  Let  I  GSj  be  the  initial  state  and  j  Conditioning  on  H,  the  sojourn  time  in  the  initial 

state,  we  have 

>  f/,2.z(c/i2+)  =  y  |//=A,Z(0)  =  |) 

J-**  ^  if  A  <  x/r,  and  Z(A  )=y  ^5,2 

=  e-'*A/,2,»,y(«  ,*).  if  A  <  x/r,  and  Z(A)=AG52 

e-**[A/i'*.y(«,*-r.  A)+  E  A^u,»(»  .»-nA  )A/,2.*.y(s  ,x)],  if  A  <  x /r,  and  Z(A  )=A  g5,. 

w  CSj 

Unconditioning  on  H ,  yields 


»A,  *A. 

A/i2,i  y(s  ,x )  =  /  e~’^dQij{k)+  E  /  «~**A/,2,*.y(s>*)^<?i»(A  ) 

0  *€5j  0 

»A. 

+  E  1  /  «'’‘Afi'»,y(«.*-fiA)dQ.*(A)  +  E  «"*A/u.«(«,*-r,A)A/,2, x)</Qit(A)l. 

*€S,  0  '  weSj 


Rearranging  yields  equation  (5). 


I 


I 

I 


I 


i 


I 

I 


i 


b)  Let  I  GSi  be  the  initial  atate  and  j  13.  Then 

E(e"'’'^,T(t)>U,i,Z{V,i+)~=j  I  Z(0)=i) 

=  E{c-*'’\T(z)>UuZ{U,+)^j  1  Z(0)=0 

+  E  E{e'*'’';T{x)>Ui,Z{Ui+)=k  \  Z(0)=i)  E{t'*"^-,Tix)>Ui2,Z{Ua+)=j  U(0)=*). 

tCSg 

Rewriting  yields  equation  (6).  Q.E.D. 

The  stochastic  properties  of  the  job  completion  process  while  the  structure-state  process  is  in  Sjj  are 
completely  described  by  the  quantities  A/i2,,(«,x)  and  A/i2,,-,/(«  ,x ),  »  i  Knowing  these 

quantities  we  can  determine  the  distribution  of  T ,  the  job  completion  time,  as  will  be  shown  in  the  next 
section. 

8.3.  The  Anal$$u  of  the  Completion  Time  th  the  pro-pri~prd  Sett  of  States 

In  this  section  we  describe  the  method  of  computing  the  LSTa  F,-  (s )  ,  i  G^i  U  U  5j.  The 
m«m  result  is  given  in  theorem  5.  First,  let  us  make  the  following  definitions.  For  i  G5, 

f.  (i-«.  (*  A. 

0 

00  •/^ 

+  E  /A/rz*(«.»)(  /  e-'*dQa(A))rf<?A) 

*€5,  0  0 

00 

+  E  /  /  e‘*‘A/,2.*(s.x-r.  A)«l<?i*(A)dG(*). 

/1C5,  0  0 

For  i  ,i  G5j 

00 

*./(•)-/  I  €-^iQii[h)dG(x) 

0  0 

00  'A 

teSfO  0 


teSjO  0 


For  iGSt  U  5, 


Miii.i(* )  *  /Af  i‘2..(«  ,*  )<*G  (* ) 
0 


and 


Wi2...y(» )  =  /A^i2.i.y(*  •*  (*  )•  J  ^ U  Sj. 

0 


Theorem  5. 

a)  The  LST»  of  job  completion  time,  F,-  (« ),  for  •  £5s  ,  satisfy  the  foUowisg  equations 

E  *.y(OA’/(*).  ‘€5,. 

i  €5g 

b)  The  LSTe  of  job  completion  time,  Ff  (s  ),  for  «  C5|  U  5]  ,  are  given  by 

Ff  (• )  =  A/ia,f(« )  +  E  (•)■  *€S|U52. 

y€5. 

Proof. 

a)  Let  t  C5g  be  the  initial  state.  Conditioning  on  H ,  the  sojourn  time  in  the  initial  state,  and 
the  initial  work  requirement  B ,  we  have; 

£(«-*■  lB=x,//-fc,Z(0)=i) 

e~"^^'  ,  if  A  >  */r, 

e-^  E  A'*'  («  )  -t-  «■**  E  (A^ «,*(«  >*  )  +  E  Af i'2,*.y(«  ,*  )F/  {»  )] 

*£J,  *€Sj  y€5, 

+  «■**  E  IA/i2.*(*  >*-»■•  *)  +  E  A/,2,*.,(<  ,v-r,  A  )Fy  (s)],  if  A  <  x /r, 

*£5,  j£5. 

Unconditioning  on  H ,  yields 

E(«-^  |B-=x,Z(0)>=i) 

*=c-“/'*(l-(?,(x/r,)) 


+  i;  /  e-**  /•*  (.V<?*(A) 

tes,  0 

•Ai 

+  £  /  *”**  (^ia,*(*  .»)  +  £  (*)]^^?<*(*) 

*€5,  0  ye5, 

•A. 

*€5,  0  yes. 

Unconditioning  on  the  initial  work  requirement,  B ,  and  rearranging  yield  equation  (7). 

b)  Let  I  €5i2  be  the  initial  state.  Conditioning  on  the  initial  work  requirement,  B ,  we  have 

E(e-^  \B=t,Z{0)^i) 

^E(e-^  ;T  <  l/,2  |  ,7(0)=i ) 

+  j:  >  Uu;Z(c;,2+)=y  |B*x.z(o)-.)£(e-*- 1  z(oHy) 

yes, 

=  A/,‘a.f(s  ,* )  +  £  ,t )  F/  («  ). 

yes. 

Unconditioning  on  B ,  yields  equation  (8).  Q.E.D. 

Theorem  5  provides  the  LST  of  the  job  completion  time  given  the  initial  state.  We  recall  that  the 
job  starts  being  serviced  upon  the  transition  of  the  structure-state  process  to  the  initial  state. 

The  procedure  to  compute  the  distribution  of  the  job  completion  time  can  now  be  stated  as  follows: 

1.  Obtain  A/i,/1[s  ,w ),  •  C5|  by  solving  equations  (1). 

2.  Invert  A/i_,‘(s  ,u> )  with  respect  to  w  to  obtain  Mi,,-(s  ,s ),  i  GSj. 

3.  Obtain  Mi,,*y(s  ,w ),  •  C5|  ,  j  ^  5]  by  solving  equations  (2). 

4.  Invert  ,u> )  with  respect  to  w  to  obtain  ,* ),  i  G5  j,  j  ^  Sj. 

5.  Obtain  Afi2,{(s  ,x  ),  i  CSj,  by  solving  equations  (3). 

6.  Compute  A/|2,,(s  ,x ),  i  C5],  by  using  equations  (4). 

7.  Obtain  ,x ),  i  €.Sq,  j  ^  S^,  by  solving  equations  (5). 


8.  Compute  ,x ),  i  GSi,  j  GS^,  by  using  equations  (6). 

9.  Obtain  Fi  {» ),  t  C^s,  by  solving  equations  (7). 

10.  Compute  F{  (s ),  i  C5|  U  S^,  by  using  equations  (8). 

11.  Compute  F  (s)» 

1-1 

12.  Invert  F  (s  )  to  obtain  F (I  )*=P(T  <  I ). 

It  should  be  noted  that  for  a  general  SAfP  the  problem  of  inversion  is  quite  difficult.  In  the  next 
section  we  illustrate  the  procedure  for  a  simple  system  with  a  structure^tate  process  that  b  a  CTMC .  In 
section  5  we  also  consider  an  application  with  a  semi-Markovian  structure-state  process. 


4.  A  Sfeeist  Cue:  A  hierkamam  Struetare-Suu  Praeeee 

Consider  a  astern  with  four  structure  states  {0,1>2,3}.  State  “O'*  b  an  absorbing  failure  state. 

Sute  “1“  b  pre ,  state  “2“  b  pri  and  state  “3“  b  pri.  Thus,  5,  *-  {1},  5,  ->>  {2},  5,  -  {3}.  The 

structure-state  process  b  a  CTMC  with  being  the  rate  of  transition  from  state  *i*  to  state  */* 

(y  y4  I ).  The  reward  rates  in  the  states  “1”,  “2”  and  “3”  are  fj.rj  and  rj,  rwpectively.  As  state  “0” 

s 

b  absorbing,  we  have  Xq/ «0,  for  all  j ,  and  vq  — =  0.  Let  X,-  ^  X,-,-  be  the  total  transition  rate  out 

f-o. 

>  ^  • 

of  state  *  I " . 

We  follow  the  procedure  stated  in  section  3  to  solve  thb  example. 

Step  1  8  i.  Obtain  Mi/{$  ,w  )  from  equation  (1). 

77^U-Ci’  (*+»‘i»)l 

$  |W 

_  V-  u _ h _ 1  _ 1; — . 

s-ffiw‘  s-t-rlW■^X^  s-^rIW+Xl 


Inverting  with  respect  to  w  ,  we  get 


where  we  have  deHned 

e,  -  .  .•  -  1,2,3. 

Step  S&l  Obtain  A/u's  («  ,w )  and  A/,*i*8  (e  ,w )  from  equation  (2). 


Inverting  with  respect  to  w ,  we  get 


Stepe  5  &  t.  Obtain  Af  I'j,!  («  ,t )  by  aolving  equatkm  (3). 

^ lifi  (•  >* )  “^  J s  (•  ,* )  ■^  ^aa  (•  »*  185  (•  »* ) 

#a  (*.*) 


11-Ajb(s  ,* )] 


Now, 


»A» 

fa  (s  ,a )  =  /va))  +  /  « "**  A/i,i  (*  ,* -r jA  )d(J 2i(A  ) 

0 


■=  e  2  +  X2ir  i(e  i-<  a)/*" 


where 


r 


r,(«  +X2)-r2(*  +^i) 


*«(•>*)*  / 


(•  +^l)  0 

(,  +x,X«  +\a)  >-’■»(* 2)1- 


Thus, 


ll-Aa(4  ,x)] 


1(«  +Xi)(*  +X2)-XjjX2i]+Xi2Xjj(-^(r,(«  +X2)e  i-fa(*  +Xi)e2)J 


(«  +X,X*+X2) 


and 


w-  /  .  93  (••*) 

« iw  (*  f* )  =  .  t  , — 7 


_  {(*•  -Xaif  i)ga-f  X2if  |«  |)(«  d-XiX*  -f  X2) _ 

71(7+xJxr+x5-x^5^»il4Ti^nrJ(r+x5«T^'^5*'+xI5^ 

From  equation  (4),  we  have 

Af  ia.i  («,*)*=  A/ 1,1  (•,»)  +  A/ 1,1,8  (•  ,t  )M  12,8  (•  f* ) 

“  « 1  +  ’*  ^• 

Steps  7  8  8.  Obtain  Mi2as  (*  >* )  solving  equation  (5). 

Af  18A8  (••*)“  •as(*  F* )  +  *8s(*  F*  )A/ia,8.s  (•  f*^ ) 

^  •a{«  ,*) 

“  11-»88(«f»)] 

Now, 

•A*  •/*’2 

«as(sF*)“  /  «'**^^a8(*)+  /  «“**A/i,i,8(*F*-v»A)d98i(A) 

0  0 

•A*  V  \  '^'s 


0  0 

0 

«o  •/»•» 

+  ^33/  /  ^  (•  I* dG{s) 

0  0 


+  X„/  /  rfG(x) 

0  0 

““  **  +  /  '  / (1“*  s)-^  12^  (<  I*  (* )  +  ^si/  /  A/i2_j  («  ,x-rsA  )rfA  iG{z) 

I*  +AjJ  {,  5  0 

with 


X.'  =  C‘((*+X,.)x/r.)  =  7x‘<'''"*’'^'‘rfG(x),  .•  =  1^.S 


and 


00  00  •/'• 

*«(•)»/  /  e-^A/«A»(«,*W««(A)^0(*)  +  /  /  ;«(•,, -r,A)dQ„(A)dC(x) 


0  0 


0  0 


00  •  A»  00  •/'• 

Aa/  /  *"^  A^isa* (•  »* 4G{z)  +  \$iJ  J  ^ A/|2,i,j{# ,x-rjfc )dA  dG(x) 
00  00 


o»Ao 


7;^/(l-e,)A/,'a,2.,(3,x)dG(x)  +  X„/  /  Af,a.u(.  ,x-r,A )«  dG(*). 

V*  0  0  0 


Thus, 


(•) 


fa(«) 

ll-A»(«)] 


00  00*  A* 

(«+Xj)ej  +X32j(l-e  j)A/j2,,(«  ,x  )dG(x  )+X8i(* +Xj))J  /  *"^  Afia,i  (•  jX-rjA  )dA  dG(x) 
_ 0 _ o_o _ 

00  00*  A» 

(•+Xj)-XjjJ(l-ej)Jl/ia^aj(#,x)dG(x)-Xj,(*+Xj)/  /  *"*  ^  (•  ,x -r  jA  )dA  dG(x) 

0  00 


From  equations  (8),  we  have 


ft  “  A/ia,i  (s )  +  A/ 18,14  (s  )F  8  (s ) 


Fi  {•)  ^  -i-  M (•  )/*»  (• ). 


Cleariy, 

•)^’/(«) 

•  -I 


5.  An  AppKeatiou:  A  Sf$Um  wiA  Two  T^pet  of  Breakdowm  amd  Repair 

Consider  a  job  to  be  executed  on  a  computer  system.  The  system  is  subject  to  three  types  of  Pois¬ 
son  fmlures: 

1)  Preemptive-resume  (pre )  failures,  which  occur  only  when  the  system  is  operating  normally,  i.e.,  in 
state  “1”,  at  a  rate  X13.  Following  a  prs  failure  the  system  undergoes  a  breakdown  period,  of  a  gen¬ 
eral  distribution  Rj^a)  and  LST  R^  (s  )>  after  which  the  interrupted  job  resumes  service. 

2)  Preemptive-ref'^atrdifferent  (pri )  failures,  which  occur  when  the  system  is  either  operating  normally 
or  undergobg  a  pre  breakdown  period,  i.e.,  either  in  state  “1”  or  in  state  “2”,  at  rates  Xj3  and  Xjs, 
respectively.  Following  a  pri  failure  the  system  undergoes  a  breakdown  period,  ot  a  general  distri¬ 
bution  R  3(s )  and  LST  R^  (<  ),  after  which  a  new  independent  job  with  work  requirement  distribu¬ 
tion  G  (r )  is  restarted. 

3)  System  failures,  which  occur  in  any  state  of  the  system  operation,  i.e.,  states  "1”,  “2”  and  “3”,  at 
rotes  Xio,  Xgo  and  X30,  respiectively.  These  are  absorbing  failures.  A  state  transition  diagram 
representing  the  system  is  shown  in  figure  1.  In  this  case  we  partition  the  states  as  follows  (in  accor¬ 
dance  with  section  3):  the  pre  set,  5, ={1,2},  and  the  pri  set,  S3={3}.  Note  that  Sj  is  an  empty 
set,  since  there  are  no  pri  states.  Obviously,  we  set  r|=l  and  r2=r3=0. 

In  order  to  obtain  the  LST  of  the  job  completion  time,  F\  (s  ),  we  follow  the  procedure  in  section  3. 

Let  H  be  the  holding  time  in  the  initial  state,  then 

<  s;Z(//-f)-y  |Z(0)-1),  , -0,2,3 

and 


<?,(.)- P(lf  <.  I  Z(0)  -  1)  -  +  «t2(*)  +  ««(*)• 

Let  X|  X}o  +  ^13  +  )^u>  PotmMi  prc^rty  of  the  different  types  of  failures,  it  foUows  thst 


and 


Qi  {•) 


^1 


e  +Xj 


Steps  1  &  t.  Obtun  Mi,i*  (t  ,iv )  and  M  1,2*  (*  )  from  equations  (1). 

1-Qi  (s+w) 


Af  1.1' (*+“’) 


s  +tl) 


+  <3i2(*+«’)Afu'(»,«')  , 


A^iV(*.w)*=  Qii  {tWi/i*  .•)■ 

It  fdlows  that 

Af  i‘i*(«  .»)—(•  +•  +X,-XjaQ,‘,  (s  )]*». 

Inverting  with  respect  to  « ,  we  have 

Af,;,  (s ,,)« 

and 

Af.,i  (•  .X )  =  (•  )A/m  (s  .X )  =  Q„  {s . 

Steps  S  &  4-  Obtain  Mi,\%  (s  ,w )  and  Af  1,2,3  (s  ,v )  from  equations  (2). 

Af  U,3  (•.»)*  ^  12  (x  +«'  )Af  IAS  (• .»  ), 

Af  iA*s  (•  •* )  ■“  +  ^21  (*  )Af  1,14  (•  #» )• 

V 

It  fdlows  that 


A^u.*s(*  •*) 


A|8^-At2^as  (*  ) 

wls+w+Xi-Xia^ai  (•)! 


We  also  have 


/  e-  (.+X«). 

0 


<22,  (* )  =  P (^  <  X  ■,ZiH+)  =  y  I  Z(0)  =  2).  y=0,1.3. 


Then 


(?2.(*)  rfi? 2(0 


and 


<22,  (*)=/«-“  rf«2l(») 

0 

=  /  e—  «2(*)  =  Pa‘  (•  +Xao+Xa»). 

0 

Similarly, 

<?«(*)  e-^{l-P,(01* 

0 

and 


00 

9a  (•)  =  /«-  rf9a(^) 
0 

0 


Xgs 

(s  +X20+X23) 


(1-P2 


(«  +X20+X2s)|. 


Let  X2  =  X20  +  X23  and  X3  =  Xso-  Then  substituting  in  F,  (»  ),  we  finally  get 

^  .  G  {$  +Xi-Xj2F2  (s  +X2)) 

(• +X8)[X,j(s +X2)+Xi2X23(1-F2  (• +X2))](1-C  (s+Xj-Xj2F2  (* +X2))] 

(s  +X2)[s  +Xj-X,2/?  2  (•  +X2)] 


Id  the  case  where  ^rstem  failures  may  not  occur,  i.e.,  X|o  Xx  Xjo  =  0,  and  X]g  «  X23,  the  above  sys¬ 
tem  corresponds  to  a  model  considered  by  Castillo  and  Siewiorek  |2].  In  this  case  F,  (s  )  reduces  to 


Fi  (•) 


_ (*+X|3)C  (« -fXi-Xiaflg  («+Xia)) _ 

(• +Xij)-Xij/?3  («  )[1-G  (j+Xi-Xi2i?2  (* +Xis))] 


Furtherm(x«,  if  the  job  service  requirement  is  deterministic  and  equal  to  x  (i.e.,  G  (*)  =  «"*• )  then 
Fi  (s  )  - - - - — - - - . 

(•+X,s>-X,sfls 

The  mean  completion  time  of  a  job,  F  ( T ),  is  given  by 

E(T)  =  -Fi'  (0)  =  I— -t-  E(Fs)] 

Xj3 

where  F(i?s)  =  -Us  (0),  is  the  mean  of  the  breakdown  period  after  prd  failures. 


6.  Comeltuiamt 

We  have  presented  a  general  model  for  the  analysis  of  job  completion  time  on  a  system  subject  to 
changes  in  its  structure  due  to  different  events.  The  system  behaviour  is  described  by  a  8emi>Maikov  pro¬ 
cess.  A  change  in  the  system  operation  preempts  job  service  which  may  later  be  resumed  or  restarted 
(with  or  without  resampling)  at,  possibly,  a  different  service  rate.  We  have  derived  a  procedure  to  obtain 
the  distribution  of  the  job  completion  time.  A  closed  form  expression  for  the  Laplace  Stieltjes  Transform 
of  the  job  completion  time  is  obtained  in  special  cases.  In  the  general  case  we  resort  to  numerical  tech¬ 
niques;  this  is  an  open  problem  for  further  research.  In  this  paper  we  have  restricted  our  attention  to  the 
execution  of  a  single  job  on  the  system.  The  obvious  extension  to  the  case  where  an  additional  delay  may 
be  experienced  due  to  queueing  is  being  investigated. 
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